Current-induced torques due to compensated antiferromagnets 
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We analyse the influence of current induced torques on the magnetization configuration of a fer- 
romagnet in a circuit containing a compensated antiferromagnet. We argue that these torques are 
generically non-zero and support this conclusion with a microscopic NEGF calculation for a circuit 
containing antiferromagnetic NiMn and ferromagnetic Co layers. Because of symmetry dictated 
differences in the form of the current-induced torque, the phase diagram which expresses the depen- 
dence of ferromagnet configuration on current and external magnetic field differs qualitatively from 
its ferromagnet-only counterpart. 
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Introduction — Current-induced torques in noncollinear 
ferromagnetic metal circuits were predicted over 10 years 
ago [E[2]> and have since been the subject of an extensive 
and quite successful body of experimental and theoretical 
research. Almost all studies of current-induced torques 
consider either their role in ferromagnetic (F) spin valve 
] or their influence on magnetic do- 
J S El ELI EDI- In both cases, the 
current-induced torques can be understood as following 
from the transfer of conserved spin angular momentum 
from current-carrying quasiparticles to the magnetic con- 
densate, hence the term spin-transfer torque. It has re- 
cently been predicted that current-induced torques are 
generically present whenever non-equilibrium quasiparti- 
cles interact with non-collinear magnetic order parame- 
ters, even [l3l[l3 | in circuits containing only antiferromag- 
netic (AF) elements. Experiments [15| have established 
a dependence of unidirectional exchange bias fields on 
current, providing indirect evidence that current-induced 
torques are present in AFs. In this Letter we analyse the 
influence of current-induced torques on a F thin film in 
a circuit containing a compensated AF. Because of a key 
difference in symmetry compared to purely F spin-valve 
circuits, we find that the phase diagram which expresses 
the dependence of the magnetic configuration of the F on 
current and external magnetic field differs qualitatively 
from the familiar F only spin- valve phase diagram [16j. 
In particular, we find that transport currents can drive 
the F to a stable steady state with magnetization per- 
pendicular to the AF layer moments. In the following 
paragraphs we argue on symmetry grounds for the form 
of the current-induced torque, and explore its robustness 
by performing a fully microscopic current-induced torque 
calculation for a circuit containing Co and NiMn layers. 
We then turn our attention to the construction of the 
F state phase diagram implied by equations of motion 
which include the current-induced torque term, and con- 
clude with a discussion of experimental implications. 
Current-induced Torques due to Compensated 
Antiferromagnets — The total current-induced torque 
acting on a F nanoparticle can always [l| be expressed in 
terms of the difference between incoming and outgoing 
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FIG. 1: Current-induced torques due to a compensated an- 
tiferromagnet. The arrows above the structure indicate the 
electron flux spin direction. The white arrows indicate the 
ensuing current-induced torques on the FM. 



spin currents. The presence of a ferromagnet will in 
general induce a nonzero spin current at the AF-F 
interface. When spin-polarized electron flux from an AF 
with orientation uaf enters a F with orientation hp, 
the spin current entering F will have some component in 
the fiAF direction. It follows that, just as in the familiar 
case where both materials are F, a current-induced 
torque will act in the plane defined by uaf and fiF, as 
illustrated in Fig. ([1]). (Out of plane torques are also 
non-zero but tend to be much smaller.) Spin-invariance 
of the overall circuit implies that the in-plane torque 
must be an odd function of the angle 9 between hp and 
fiAFi and that it can therefore be expanded in terms of 
a sin-only Fourier series, vanishing for both parallel and 
antiparallel collinear configurations. Most AF materials 
used in magnetoelectronics are fully compensated, i.e. 
the spin-density sums to zero (or nearly so) in every 
lattice plane perpendicular to the current direction. 
In this case, reversal of the AF moment direction is 
equivalent to a lateral translation which cannot influence 
the current-induced torque. It follows that in the 
compensated AF case the torque is invariant under 
9 —> 9 + 7r, restricting its Fourier expansion to terms 
proportional to sin(2n#). The torque therefore vanishes 
when hp is perpendicular to uaf, and undergoes a 
sign change for 9 — > it — 9, as illustrated in Fig. |T]). 
The property that the torque acting on a F due to a 
compensated AF vanishes not only for collinear but also 
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for perpendicular orientations is primarily responsible 
for the novel current-induced torque phase diagram that 
we discuss below. 

Current- Induced Torques for Co/NiMn — We employ a 
non-equilibrium Green's function (NEGF) approach 17 1 
for microscopic calculations of magneto-transport prop- 
erties and current-induced torques. Quasiparticle Hamil- 
tonians are constructed using density functional theory 
within the local spin density approximation (extended to 
allow noncollinear spin configurations), norm-conserving 
pseudopotentials, and an s, p, d single-zeta atomic or- 
bital basis set. The induced torque per-current can be 
calculated atom by atom 17j : 
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The right-hand-side of Eq. |T]) expresses the misalign- 
ment between the non-equilibrium spin-density, rn 1 ^ p, 
and the spin-dependent part of the exchange-correlation 
potential, A a ^. Here T{ep) is the transmission probabil- 
ity, fc|i labels transverse channels, a, (3 are orbital labels, 
and a is summed only over orbitals centered on the atom 
of interest. When the many-body Hamiltonian is spin ro- 
tationally invariant, the current-induced torque on each 
atom is equal to the net spin flux out of the atom. 

We apply this approach to a system with a single in- 
terface between antiferromagnetic NiMn and ferromag- 
netic Co. The crystal structure of NiMn is face centered 
tetragonal, with Ni and Mn layers alternating in the (001) 
direction 18[ . The Ni atoms are approximately nonmag- 
netic, while the Mn atoms form a compensated antifer- 
romagnetic 2-dimensional lattice within each plane (See 
Fig. 101). In our calculation, we use a = 3.697 A, with 
a c/a ratio of 0.9573 for NiMn and, following Ref. [l9j ]. 
a lattice-matched tetragonal structure for Co with a c/a 
ratio chosen to conserve its experimental atomic volume. 
The results shown here are for current in the (001) di- 
rection, perpendicular to the interface between Co and 
Ni terminated NiMn. The current through the interface 
has a polarization P = (T T -T x ) / '(T T +T l ) = 6.4% when 
the F and AF moments are collinear, with the larger con- 
ductance for the ferromagnet majority spins. To evaluate 
the current-induced torques present in the system, we ro- 
tate the Co layer magnetization orientation by an angle 
9 with respect to the NiMn moment direction and use 
Eq. ©. 

Fig. ([!]) shows the total torques acting on the AF 
and F order parameters as a function of 9. The torque 
acting on the F closely follows the form anticipated on 
symmetry grounds above. We associate the small torque 
at 9 = 7r/2 with weak ferromagnetism which is induced 
in the top (Ni) layer of NiMn. In NiMn only the differ- 
ence between the torques on the two sublattice of the AF 
drives the order parameter. The current-induced torque 
tends to drive the orientation of downstream material 
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FIG. 2: Illustration of the NiMn-Co interface model. 



(AF or F) parallel with that of the upstream, and to 
drive the upstream material orientation perpendicular to 
the downstream (so that for electron flow from AF to F, 
the F tends to align to AF, and the AF tends to become 
perpendicular to F, within their common plane). 

We have also considered Mn terminated NiMn adja- 
cent to Co. In this case the last Mn layer acquires a net 
magnetic moment in the direction of Co. The current- 
induced torques do not show as clean of a sin 29 behavior, 
but a combination of sin9 and sin29. We conclude that 
the absence of odd n sin n9 torques is closely tied to the 
degree of compensation at the AF interface. 



STT/I acting on AFM and FM 



■J) 




-0.25 



-0.5. 



45 90 135 180 
Relative orientation 



FIG. 3: Current-induced torques per current acting on the or- 
der parameters of the AF and F layers vs relative orientation. 
In the AF the order parameter is driven by the differences 
between torques on opposite sublattices. Units are ps/e. 



Phase Diagram for a pinned antiferromagnet — We now 
consider the implications of this new form of current- 
induced torque for systems with the usual thin film geom- 
etry, assumming for the sake of definiteness that the AF 
moment direction haf is pinned and lies in the plane, and 
that the external magnetic field H is applied in the same 
direction. We use a spherical coordinate system for the F 
moment direction, taking uaf as the polar direction and 
the x direction as the film normal (the demagnetizing 
field is denoted by Hj). We assume that a non- magnetic 
spacer layer is placed between the F and AF layers so 
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that exchange interactions negligible. Just as in the pure 
F case, a spacer layer is not expected to have a large im- 
pact on current-induced torques. We also omit easy-axis 
anisotropy; its inclusion wouldn't substantially change 
the picture described below. With these ingredients the 
polar and azimuthal torques acting on the ferromagnet 
are: 



r 
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Here we have parameterized the current-induced torque 
by Hci and chosen a sign convention in which Hci < 
when it favors perpendicular alignment. Steady-state so- 
lutions satisfy = Tg = and are stable for small devi- 
ations when Gilbert damping is included. We have deter- 
mined the stability regions of the steady state solutions 
discussed below by following the procedure described in 
Ref. 0. We present all of our results in terms of the 
dimensionless fields h = H/Hd and hci — Hci/Hd- 

In the absence of the current induced torques, the mag- 
netization simply lines up with the magnetic field applied 
in the easy plane. The influence of the current-induced 
torque on F is particularly dramatic for Hci < 0. Be- 
cause the torques then tend to push the magnetization 
perpendicular to haf, the field-aligned solution is stable 
only for: 
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where a is the Gilbert damping parameter. For suf- 
ficiently strong currents and weak external fields, a 
perpendicular-to-plane steady state becomes stable: 



= —2hcih + Tin 



(4) 



These equations have been derived assuming that h and 
hci are small. In the above n is an even integer for 
solutions which point approximately in the +x direction, 
and an odd integer for the —x direction. The region of 
stability for this solution is: 



hci < — 2 



h sin h — 2 cos 2 h 
hsmh — cos2/i 



(5) 



where the the fraction on the r.h.s. of the above in- 
equality must be negative, implying that \h\ < 0.608, or 
equivalent^ \H\ < /j, M s (0.608). 

The stability of this counter-intuitive stable steady 
state is explained in Fig. Q. This figure illustrates 
the situation when the excursions from the easy plane 
are small. For simplicity we first consider no external 
field. In the absence of the current-induced torque a 



small fluctuation out of the easy plane would initiate pre- 
cession about the hard axis which damps back into the 
easy plane. The presence of the sin 29 torque, however, 
drives the magnetization rh perpendicular to uaf within 
their common plane. As rh precesses around the hard- 
axis, this torque vector has a component which points 
out of the easy plane. If the angle between Haf and the 
in-plane component of m is /3, the magnitude varies as 
T x = 2Hcim x sin 2 (3, as shown in the figure. The crucial 
point is that this torque is always positive throughout 
the precession. When this torque exceeds the damping, 
the out-of-plane configuration is stabilized. The eventual 
out-of-plane orientation can be +x or — x depending on 
the direction of the initial fluctuation out of plane. The 
presence of an applied field changes the trajectory of the 
magnetization upon excursions from the easy-plane. For 
a sufficiently large applied field, the torque is unable to 
stabilize the out-of-plane configuration, and no steady 
state is reached. 




FIG. 4: Saddle shape illustrates the out-of-plane torque vs j3 
for small excursions of the magnetization orientation rh from 
the easy plane. The out-of-plane torque is always positive. 



Interesting new steady states can in principle also be 
induced by the current-induced torque for Hci > 0. For 
\H\ < Hd, the steady state stability analysis identified 
configurations in which the magnetization is approxi- 
mately anti-aligned with the applied field: 



= cos 



(-h) 



-2h C ih. 



(6) 



which is stable for the range of applied fields and currents 

2-h 2 



l CI > 77 
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For Hci > and \H\ > Hd, the equilibrium solutions 
are m z = ±1. In this case the stability condition for the 
magnetization anti-aligned with the field is: 



i^ci 



> 



h\ 



(8) 



These anti-aligned states occur only if the magnetization 
is initially nearly anti-aligned with the applied field. The 
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reason for their stability is that this form of the current- 
induced torque does not distinguish between +z and — z 
- it merely tends to make to direct the F to the near- 
est available z-axis, even if it's opposite to the applied 
field. The region for such a solution is shown in Fig. 
labelled ±z. This misaligned steady state may not be 
experimentally relevant however because it occurs only 
when the magnetization is initially nearly anti-aligned to 
an applied field of finite magnitude \H\ > i?d\/l/3. 

Fig. ([5|) shows the x and z components of the magne- 
tization as a function of applied field and current, deter- 
mined numerically. We have taken the damping a = .01. 
Also shown is the magnitude of the power spectrum peak 
of z{t) (labelled "Pz") - a nonzero value indicates a pro- 
cessing solution. Also shown is the stability boundaries 
defined by Eqs. ([21 EJ El [5]) . The numerics verify the sta- 
bility of the unusual out-of-plane and field-anti-aligned 
solutions. The conversion of the dimensionless hci into 
a real current density for a material with demagnetiza- 
tion field of 1 T is J = (hcit) x 3.8 • 10 9 A/cm 2 , where t 
is the thickness of the F layer in ran. 




Stability Boundaries 
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FIG. 5: Magnetic configuration (M X ,M Z ) and peak of power 
spectrum P z (arbitrary units) versus applied field and cur- 
rent. Also shown is stability boundaries found analytically 
(the labels ±x, ±z refer also to solutions which point approx- 
imately in these directions). The stability boundary plot also 
shows the reduced out-of-plane solution space for negative to 
positive field sweep with a dashed line. 



The data for each (h, hci) point of Fig. ([5]) is obtained 
beginning from an initial condition close to the solution 
given by Eq. ([4]). These equilibrium solutions are not 
universal attractors, and are attained for a subset of ini- 
tial conditions. To see the effect of initial conditions, we 
have also swept the applied field from negative to posi- 
tive for each applied current, using the slightly perturbed 
final coordinates of a trajectory as the initial condition 



for the next value of applied field. The out-of-plane so- 
lution space is reduced, shown by the dotted line in Fig. 
([5]) in the stability boundaries plot. 

We now comment on the experimental possibilities of 
seeing these effects. In the preceding analysis, we assume 
that the AF is fixed. This can be accomplished by placing 
a large F adjacent to the AF, so that the AF is pinned via 
the exchange bias effect (the overall stack structure would 
be pinning F - AF - spacer - free F) . The presence of this 
pinning F may influence the dynamics of the free F, but 
its signature should be very distinct from the influence of 
the AF layer on the free F. The orientation of the free F 
should be observable from magnetoresistance effects with 
the pinning ferromagnet. 

A virtue of the out-of-plane F configuration is that the 
surface of the AF need not be single domain for its ob- 
servation. As long as the magnetization of the AF is 
compensated and points in t he p lane (which is the pre- 
ferred direction for NiMn 2l|, |22|), different orientations 
of domains at the AF surface should cooperatively push 
the F out of the plane. The encouraging aspect of this 
proposal is that the signature of the AF current-induced 
torque is so unique, helping to provide a distinguished 
characteristic for its observation. 

Acknowledgments — We would like to acknowledge 
very helpful conversations with Maxim Tsoi and OUe 
Heinonen. This work was supported in part by Seagate 
Corporation and by the National Science Foundation un- 
der grant DMR-0606489, and the computational work 
was supported by Texas Advanced Computational Cen- 
ter (TACC). 



[1 
[2 
[3] 
[4] 

[s: 

ir 

[8] 
[9] 

[io; 
[ii 

[12 
[13 

[14 
[15 
[16 
[17 
[18 

[19 
[20] 



Slonczewski, J. Magn. Magn. Mat. 62, 123, (1996). 

L. Berger, Phys. Rev. B 54, 9353 (1996). 

M. Tsoi et al., Phys. Rev. Lett. 81, 493(E) (1998). 

M. Tsoi et al., Nature 406, 46 (2000). 

J. A. Katine et al, Phys. Rev. Lett. 84, 3149 (2000). 

F. J. Albert et al, Phys. Rev. Lett., 89, 226802 (2002). 
M. R. Pufall, W. H. Rippard, and T. J. Silva, App. Phys. 
Lett. 83, 323 (2003). 

G. S. Beach et al, Phys. Rev. Lett. 97, 057203 (2006) 
J. Grollier et al., Appl. Phys. Lett. 83, 509 (2003). 

A. Yamaguchi et al, Phys. Rev. Lett. 92, 077205 (2004). 
M. Yamanouchi, D. Chiba, F. Matsukura, and H. Ohno, 
Nature 428, 539 (2004). 

M. Hayashi et al, Phys. Rev. Lett. 96, 197207 (2006). 
A. S. Nunez, R. A. Duine, Paul Haney, and A. H. Mac- 
Donald, Phys. Rev. B 73, 214426 (2006). 
P. 
Z. 

s. 
P. 

L. 
T 



M. Haney et al, Phys. Rev. B 75, 174428 (2007). 
Wei et al., Phys. Rev. Lett. 98, 116603 (2007). 
I. Kiselev et al, Nature 425, 380 (2003). 
M. Haney et al., Phys. Rev. B 76, 024404 (2007). 
Pal et al., J. App. Phys. 39, 538 (1968). 
C. Schulthess and W. H. Butler, J. App. Phys. 
7225 (1998). 

Ya. B. Bazaliy, B. A. Jones, and Shou-Cheng Zhang, 



83, 



Phys. Rev. B 69, 094421, (2004). (1959). 
[21] A. Sakuma, J. Mag. Magn. Mat. 187, 105 (1998). 
[22] J.S. Kasper, J.S. Kouvel, J. Phys. Chem. Solids 11, 231 



